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Abstract 



We derive a closed formula for the matrix elements of the volume operator for 
canonical Lorentzian quantum gravity in four spacetime dimensions in the contin- 
uum in a spin-network basis. 

We also display a new technique of regularization which is state dependent but 
we are forced to it in order to maintain diffeomorphism covariance and in that sense 
Q\ | it is natural. 

We arrive naturally at the expression for the volume operator as defined by 
. Ashtekar and Lewandowski up to a state independent factor. 

On 

"o 1 Introduction 

cr 



The volume functional V for Lorentzian canonical quantum gravity plays a quite impor- 
tant role in the quantization process : if one follows the approach advertized in |lj and 
chooses a representation in which traces of the holonomy along closed loops for a certain 
SU (2) A l a connection can be promoted to basic configuration operators, then it becomes 
technically very hard to define physically interesting operators such as the Euclidean and 
Lorentzian Hamiltonian constraint operators, an operator corresponding to the length of 
a curve, an operator corresponding to the ADM energy, matter Hamiltonian operators 
and so forth. The common reason is that all these operators involve only the co-triad 
e l a of the intrinsic metric of the initial data hypersurface £ polynomially which, however, 
cannot be written as a polynomial in the momentum conjugate to A l a and therefore does 
not seem to make any sense as a quantum operator. 

Fortunately, there is a trick available : as proved in 0, |3|], the total volume of £ is 
the generating functional of the co-triad meaning that one gets e l a from V essentially by 
functional derivation. Therefore, if we write Sts a e l a oc {h s , V} + o(5t 2 ) where s is an edge 
with parameter length St, h s is the holonomy along s and {., .} denotes Poisson brackets 
then it would be sufficient to quantize V in such a way that (gauge-invariant) functions 
of a finite number of holonomies form a dense domain. Such a satisfactory quantization 
of V which has precisely this property was indeed given already in the literature ||, |[]. 
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Therefore, the strategy outlined above and in 0, |3| can be implemented and, in partic- 
ular, was successfully employed to construct various operators listed above ||, ||, || [7|, || . 
However, the explicit computation of the spectrum of most of those operators was not 
performed so far. The reason is that the spectra of all these operators are clearly largely 
determined by the spectrum of the volume operator whose precise spectrum is fairly 
unknown except for the most simple situations || [TO) . 

The present paper is devoted to the derivation of the complete set of matrix elements 
of this operator in closed form from which the spectrum can be obtained straightforwardly. 

In section 3, after introducing the notation in section 2, we develop a novel technique 
to regularize the volume operator which is state dependent very much in the same way 
as the regularization displayed in the second reference of ||. This part of the paper is 
a continuation of the work on the point splitting regularizations of the volume and area 
operators as displayed in Jl|, ||. The derivation is rather short and gives a diffeomorphism 
covariant end result. We are unambiguously led to the Ashtekar-Lewandowski volume 



operator defined in [18| up to a factor of y27/8. 

In section 4 we repeat the argument given in || |9j which shows that the spectrum 
is discrete and that its computation reduces to linear algebra, in the sense that we just 
need to determine separately the eigenvalues of an infinite number of finite-dimensional, 
positive semi-definite and symmetric real-valued matrices. 

Finally, in section 5 we derive the explicit formula for an arbitrary matrix element in 
a spin-network basis and compute the eigenvalues of the associated matrices in a number 
of simple cases. It should be noted that the formulae obtained, albeit quite complicated, 
are the starting point for a numerical evaluation of the spectrum and are in particular 
useful to determine the large spin behaviour of the volume operator. 



2 Notation 

Let us first introduce the notation. We consider the Hamiltonian formulation of the Pala- 
tini action of four-dimensional Lorentzian vacuum gravity. 

Denote by q a b the intrinsic metric on the initial data hypersurface £ and let e l a be its 
co-triad where i,j,k,.. denote su{2) indices. The densitized triad Ef := det((e^))e" 
multiplied by k being Newton's constant, is the momentum conjugate to K l a : = 

sgn(det((e^)))i ; r a fee^ where K a b is the extrinsic curvature of S. Upon performing a canoni- 
cal point transformation |jnj |12[ on the phase space coordinatized by (K l a , E^/k) we arrive 



at the chart (A l a / k : = (T* + iQ)//«, E" 1 ) where T l a is the spin-connection of e\. It is easy 
to see that the configuration variable A l a is a su(2) connection. 

The virtue of casting general relativity into a connection dynamics formulation is as 
follows. 

1) First of all, the phase space of general relativity is now embedded into that of an 
su{2) gauge theory, thereby opening access to a wealth of techniques that have been 
proving useful in quantizing those theories. Probably the most important technique 
is the use of traces of the holonomy along piecewise analytic loops in S, so-called 
Wilson-loop functionals, to coordinatize the space of smooth connections modulo 



gauge transformations A/Q |T3|, |14| 
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2) Using those Wilson loop functionals one can construct [|T5| , |IE , |T7| a quantum con- 
figuration space of generalized connection modulo gauge transformations A/Q and 
find natural, a-additive, diffeomorphism invariant, regular Borel probability mea- 
sures /i thereon which equips us with a Hilbert space structure 7i := L^A/Q, dfio). 
Of course, for this to make sense we have to make the following assumtion : Wilson 
loop functionals continue to make sense as operators in the quantum theory. We 
will adopt this viewpoint in the sequel in the hope to capture one of the physically 
interesting phases of quantum general relativity. The interested reader is referred 
to [|l| and references therein for further details. 

The following paragraph is supposed to equip the reader with a working knowledge of the 
tools associated with A/Q. 

We begin by defining cylindrical functions. Since the gauge invariant information con- 
tained in a connection is captured by finite linear sums of products of Wilson loop func- 
tionals we may label gauge invariant functions by piecewise analytic graphs 7 which are 
just the union of the loops involved. Such a graph consists of maximally analytic pieces, 
called the edges e of 7 and the edges meet in the vertices v of 7. A function / on A/Q 
is said to be cylindrical with respect to a graph 7 iff it can be written as / = / 7 o p 7 
where p 7 (A) = (h ei (A), .., h en (A)) and where ei, .., e n are the edges of 7. Here h e (A) is the 
holonomy along e evaluated at A e A/Q and / 7 is a complex valued function on SU(2) n . 
Since a function cylindrical with respect to a graph 7 is automatically cylindrical with 
respect to any graph bigger than 7, a cylindrical function is actually given by a whole 
equivalence class of functions / 7 . In the sequel we will not distinguish between this equiv- 
alence class and one of its representants. 

We say that a function cylindrical with respect to a graph 7 is of class C n if it is of class 
C n with respect to the standard differentiable structure of SU{2) n and denote this class 
of functions by Cjl^(A/Q), Cy\ n (A/Q) := U jEr Cy^(A/Q) where T is the set of piecewise 
analytic graphs in E. Thus we have a differential calculus [JT^] on A/Q. 



Finally, the measure //o is the faithful projective limit ||T9| , |20| of the self-consistent pro- 
jective family (/i 0j7 ) 7g r defined by 

/ dfio(A)f(A) := / dfio t y(A) / 7 (p 7 (A) ) := / d/j lH (gx)..diJ l H(gn)f 1 (gi,-,gn) 

JA/Q J A/Q JSU(2) n 



where gi := h ei (A). This equips us with an integral calculus on A/Q. 



3 Regular izat ion of the volume operator 

In this section we present a new regularization of the volume operator which uses the 
methods derived by Abhay Ashtekar and Jurek Lewandowski in fTB| , [|, Most of the 
issues involved in the regularization as presented here are borrowed from their work. 
The advantage of our procedure is that it is rather short and compact. 

Let R C E be an open, connected region of E. We have the identity 

\e ah ^ k E«E)El = det((£f)) = det(M) = [det(«))] 2 =: det(g) > (3.1) 
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and can thus write the volume of the region R as measured by the metric q a b as follows 



V(R) := j R d 3 xy / de%) = ^ ^e^EfEpt . (3.2) 

The next step is to smear the fields Let Xa{p, x) be the characteristic function in the 
coordinate x of a cube with center p spanned by the three vectors Aj = AjfTj(A) where 
Hi is a normal vector in the frame under consideration and which has coordinate volume 
vol = A1A2A3 det(ni, H2, n$) (we assume the three normal vectors to be right oriented). 
In other words, Xa(p, x) = nf=i — \ < n^x — p > |) where < ., . > is the standard 
Euclidean inner product and 8(y) = 1 for y > and zero otherwise. 
We consider the smeared quantity 

' f j3_ f j3_. f j3_._ , /_ _\_ . /_ x + 2/n 



E(p, A, A', A") := / d 3 x / d 3 y / d 3 z X AMxA'(p, 

vol(AJvol(A'Jvol(A") js js js 2 



x 



XXA"(P, ^|±£ )^ eafa ^*^i( a .)^(y)^( g ) . (3.3) 

Notice that if we take the limits Aj, A[, A" — > in any combination and in any rate 
with respect to each other then we get back to (|IJ) evaluated at the point p. This holds 
for any choice of linearly independent normal vectors n^n'^n^. The strange argument 
(x + y + z)/3 will turn out to be very crucial in obtaining a manifestly diffeomorphism 
covariant result. We will see this in a moment. 
Then it is easy to see that the classical identity 

V(R) = lim lim lim / d 3 pJ\E(p, A, A', A") I (3.4) 
holds. We could introduce absolute values because we have the classical identity det(g) = 



I det(g) I and the modulus is necessary for the regulated quantity ( |3.3| ) is not non- negative 
anymore. 

The virtue of introducing the quantities ( |3.3j ) is that they can be promoted to quantum 
operators which have the dense domain Cyl 1 (^4/^). To see this, notice that due to the 
canonical brackets {A\(x), Ej(y)} = K5 b a 5 % j5^ (x,y) we are naturally led to represent the 

operator corresponding to Ef by Ef(x) = —i£p5/5A l a (x) where £ p = \/hn is the Planck 
length. 

Now let be given a graph 7. In order to simplify the notation, we subdivide each edge e 
with endpoints v,v' which are vertices of 7 into two segments s, s' where e = s o (s') _1 
and s has an orientation such that it is outgoing at v while s' has an orientation such that 
it is outgoing at v' . This introduces new vertices s D s' which we will call pseudo- vertices 
because they are not points of non-analyticity of the graph. Let E(j) be the set of these 
segments of 7 but V(j) the set of true (as opposed to pseudo) vertices of 7. Let us now 
evaluate the action of E^(p,A) := l/vol(A) / s d 3 xxA(p, x)E? (x) on a function / = p*/ 7 
cylindrical with respect to 7. We find (e : [0, 1] — ► E; t — > e{t) being a parametrization 
of the edge e) 

if fi 

E«(p,A)f = -f- £ / dt XA (p,e(t))e a (t)x 

vol(A) ^~^!_.^Jo 



I , . d 

l dh e (0,l) 



x -ti(h e {0,t)nh e (t,l) — /n - J / 7 . (3.5) 
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Here we have used 1) the fact that a cylindrical function is already determined by its 
values on A/Q rather than A/Q so that it makes sense to take the functional derivative, 
2) the definition of the holonomy as the path ordered exponential of J e A with the smallest 
parameter value to the left, 3) A = dx a A l a Ti/2 where SU{2) 3 Tj = —icrj, Oi being the 
usual Pauli matrices, so that [rj/2, Tj/2] = tij k T k /2 and 4) we have defined tx{hd/dg) = 
hABd/dgAB, A,B,C,.. being SU(2) indices. 

We now wish to evaluate the whole operator E(p, A, A', A") on /. It is clear that we 
obtain three types of terms, the first type comes from all three functional derivatives 
acting on / only, the second type comes from two functional derivatives acting on / 
and the remaining one acting on the trace appearing in (|3.5| ) and finally the third type 
comes from only one derivative acting on / 7 and the remaining two acting on the trace. 
Explicitly we find (we mean by 9(t, t') the theta function which is unity if < t < t' < 1, 
1/2 if t = < t' < 1,0 < t < t' = and 1/4 if t < t' e {0,1} and zero otherwise. 
Likewise 0(t,t',t") is 1 if t < t' < t", 1/2 if t = < t' < t",0 < t < t' < t" = 1, 1/4 if 
= t = t' < t", < t < t' = t" and 1/8 if t < t' < t" G {0, 1}. Finally m(t, t', t") = l/2 k 
where k — 0, 1, 2, 3 is the possible number of arguments that equal 0, 1) 

i£ 6 r 

E(p, A, A', A")/ = o o, UM VX* u^ e ^ Jk I „ dtdt ' dt " X 
8 ■ 3!vol(A)vol(A')vol(A") ^[o,!] 3 

x: m a e\t') b e''(frxA(p,e(t))xA'(p, e{t) + / {t ' ] ) x 

t e,e',e"eE(-y) 

, e(t) + e'(t') +e"(t"), , , ... , n ... , . .. d 
x XA"(p, ^ K -^)m(t,t',t")tr(h e „(0,t")T k h e „(t", l) dh , Q x 



x tr^KO^Or,/,^^!)^^)^^,^^,^,!)^^) 



+ £ e"(t) a e'(t') 6 e"(t") c XA^, e"(t)) X A'(p, e " (t) + ^ x 

e',e"GE( 7 ) ^ 
X XA"(P, U V 1 ) )[9(t,t'')tT(h e/l (0,t)T l h ell (t,t'')T k h e/l (t'', V 



dh e „(0A) 

d d 
+e(t",t)tT(h e/ iO,t")T k h el >{t",t)Tih el/ (t,l) dh , Q 1 J ]tr(/i e ,(0,t / )r^ e /(t / ,l) afe Q v 

+ e'(tye'(t') b e"(t"y XA (p, e'(t)) XA ,(p, g^l+^gl ) x 

e',e"6B(7) 

, e'(t) 4-e , (i') + e"(t")x /rt , A , , „ „ 9 



3 y v e v ' ' * ' v ' 'dMM) 

x [9(t, t')tr(h e ,(0, *frM*, t'^jh^t', 

+0(t / , t)tr(^(0, t'^heit', t)nh e ,(t, 1) ^ ^ )] 
4- £ eft) tt e ff (0 5 ^(O c XA(p,e(0)XA>(p, ^ x 

e,e"eB(7) 

x Xa»(p, W V 1 ; )[g(^,Otr(/ie>>(0,0^(< / ,O^V(< ff > l) 



a/i e »(o,i)- 
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+9(f,t%(h e ,,(0X)Tkh e >it",t%h e r,(t'^) — )]tr(fe e (0,t)r^ e (t,l) 

oh e »{0, 1) 0/i e (0, 1, 

+ £ e"(t) a e"(t') 6 e"(t /, ) C XA(p, e»(t))xA'(p, ^1±^®) x 

e"£E(7) 

x X a» (p, — ^ Y ) x 

x [0(t,t / ,t ,/ )tr(/i e //(O,t)r i /ie"(*»^)' 7 i' l e"(* )nK"{t , 1) 



<9/i e »(0,l) 



+6*(t, t ,t)tT(he»(0,t)Tih e H(t,t )T k h e „{t ,t )Tjh e n{t , 1)— — — — ) 

+0(* / , t", t)tr(/i e « (0, t')rj -,h e " {t , t")r k h e „ (f, t^h^ (t,V 



dh e "(0,l) 
d 

+9(t , £, £ )tr(/i e «(0, £ )rjh e /i(t , t)Tih e "(t, t")r k h e "(t" , 1)- 



<9/i e »(0,l) 

+6>(£", £, £')tr(/i e »(0, t")r k h e n(t" , t)rih e »(t, £ )r^A^(t , 1)- 



a/i c »(o, i) y 

+ 0(£", £', t)tr(M°> t")r k h e „{t", t'^K^t', t)nhAt, ^) Qh //(Q jr )] | A 

= : [01,2,3 + 02,31 + 012,3 + 01,23 + 0x2^7- ( 3 - 6 ) 

The fact that the integrand of the terms involved in 0^2,3,01,23,02,31,0123 vanishes if 
either of the cases < £ = £' < 1, < £' — t" < 1, < £ = £' — t" < 1 occurs is 
due to the fact that in this case in 012,3,01,23,02,31,0123 we get a trace which contains 
TuTj\, VjTty, TfaT{\ contracted with e y which vanishes (to see this recall that the functional 
derivative is 

5h e {A)/5A\{x) = \ f 1 dt[l-S {3 \e(t+),x)e(t+) a h e (0,t)T l h e (t,l) 
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+ i 5 (3) ( e ( t _) } x )e(t-) a h e (0, t)r t h e (t, 1)] 

(one sided derivatives and 6 distributions). This expression is also correct if x is an 
endpoint of e (in which case there is only one term which survives in ( |3.7|) ). which results 
in the case that we consider h ei Tjh e2 instead of h e , e = e\ o x = e\ fl t2 a point of 
analyticity, in a term involving TuTj\). 

Given a triple e, e', e" of (not necessarily distinct) edges of 7, consider the functions 

w(t ,,, :=?W±4Q±4Q. (3.7) 
This function has the interesting property that the Jacobian is given by 

3 3 det ^ 3 ^"' ^0kA & tl f ) ) = e abc e(t) a e\t') b e"(t") a (3.8) 

which is precisely the form of the factor which enters all the integrals in ( |3.6|) . This is 
why we have introduced the strange argument (x + y + z)/3. 
We now consider the limit Aj, A' i: A" — > 0. 
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Lemma 3.1 For each triple of edges e, e', e" there exists a choice of vectors Hi, n' iy n" and 
a way to guide the limit Aj, Aj, A" — > such that 

I det( f }f^''l \ ) X a(p, e) X A'(p, (e + e')/2) X A»(P, (e + e' + e")/3)6 ee , e „ (3.9) 
J[o,i] :i o(t,t',t") 

vanishes 

a) ife,e',e" do not all intersect p or 

b) det( g( t e /|") ) P = (which is a diffeomorphism invariant statement). 
Otherwise it tends to 

d(x a ) * o 

l/8sg , n(det( g ( - ( -^, /e f " / ) ))pQe,e , ,e" (p) ]lj=i A", i/ere we /iai>e denoted by O ee < e «(t,t' ,t") the trace(s) 
involved in the various terms of (\3.b\ ). 

Proof] : 

If at least one of e, e', e" does not intersect p then, if we choose Aj etc. smaller than 
some finite number Ao, (|3.9|) vanishes identically since the support of the characteristic 
functions is in a neighbourhood around p which shrinks to zero with the Aj etc. 
So let us assume that all of e,e',e" intersect p at parameter value t ,t' ,tQ (this value is 
unique because the edges are not self- intersecting). Then we can write e(t) = p + c(t — t ) 
where c is analytic and vanishes at r = t — 1 = 0. We have the case subdivision : 

Case I : det(|gj^)) to = 0. 

Case la) : All of 6(0), c'(0), c"(0) are co-linear. 

Case lb) : Two of c(0), c'(0), c"(0) are co-linear and the third is linearly independent of 
them. 

Case Ic) : No two of c(0), c'(0), c"(0) are co-linear. 
Case II : det(|g^)) to + 0. 

Notice that all vectors c(0), c'(0), c"(0) are non- vanishing by the definition of a curve. 

We consider first case I). We exclude the trivial case that all three curves lie in a 
coordinate plane or line such that the determinant already vanishes for all finite values of 
the A's. Therefore there exist linearly independent unit vectors u,v,w (not necessarily 
orthogonal) in terms of which we may express c, c', c". 
In case la) we have an expansion of the form 

c(t) = au(t + o(t 2 )) + bv(t m + o[t m+l )) + cw(t n + o(t n+1 )) 
c'(t) = au(t + o(t 2 )) + b'v(t m ' + o(t m ' +1 )) + c'w(t n ' + o(t n ' +1 )) 
c"(t) = a"u(t + o(t 2 )) + b"v(t m " + o(t m " +1 )) + c"w(t n " + o(r" +1 )) 

(3.10) 

where a, b, c, a', b', c', a", b", c" are real numbers with aa'a" ^ and at least one of the Vs 
and c's being different from zero (also not for instance b = c = b' = d = 0). Furthermore 
m,m' ,m" ,n,n' ,n" > 2. The characteristic functions have support in coordinate cubes 
spanned by the vectors nj,n^,nf. Now, since u,v,w are linearly independent we may 

1 The idea that one should adapt the regularization to each triple of edges was first communicated to 
the author by Jurek Lewandowski. Its justification relies on the fact that the classical expression does 
not depend on the way we regularize. What is new here is to introduce even more parameters than in 
to be adapted to a triple which simplifies the proofs. What is also new is that we introduced the 
(x + y + z)/3 argument which makes the proof especially clear. 
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simply choose, for instance, := := v, n" := w. It follows then and from the fact 
that < x < 1 that 

Xa(p, e)xA'(p, (e + e')/2)xA"(p, (e + e' + e")/3) 
= Xa(0, c)xa/(0, (c + c')/2)xa»(0, (c + c' + c")/3) 

< 9 Al (<c,u>)6 A/i (<(c + c')/2,v>)9 A ,,(<(c + c' + c")/3,w>) . (3.11) 



From the explicit expansions of c, c', c" we conclude that ( 3.11|) has the bound 



e^it^s^it')^'!^") (3-12) 
for some sufficently large numbers 5i,5[,S". On the other hand we also see from the 

Q(x a ) 

explicit expansion of | det( g( t «fffl )| around t that it is bounded by M(\t\ k + \t'\ k + \t"\ k ) 
where M is a positive number and where k = min(m + n',m + n", m! + n,m' + n", m" + 
n,m" + ri) - 2 > 2. 

The prescription of how to guide the limit in case la) is then to synchronize Ai = A[ = 
A" = A and to take the limit A — > first. The integral is at least of order A 5 while we 
divide only by an order of A 3 so that the result vanishes. 

In case lb) we have an expansion of the form (let w.l.g. c, d have co-linear tangents) 

c(t) = au(t + o(t 2 )) + bv(t m + o(t m+1 )) + cw(t n + o(t n+1 )) 
c'(t) = au(t + o(t 2 )) + b'v(t m ' + o(t m ' +1 )) + c'w(t n ' + o(t n ' +1 )) 
c"(t) = a"v(t + o(t 2 )) + b"u(t m " + o(t m " +1 )) + c"w(t n " + o(t n " +1 )). 

(3.13) 

We now argue as above and find that the product of the characteristic functions can be 
estimated by 

05 1 A 1 Ai (t )0qA» (t ) 

while the determinant can be estimated as above just that k is now given by k = 



min(m, m', m", n, n', n") — 1 > 1. 

The prescription is now Ai = A[ = A'^ =: A — > first and we conclude that the integral 
is at least of order A 4 while we divide again only by A 3 such that the limit vanishes. 
In case Ic) finally we have an expansion of the form 

c(t) = au(t + o(t 2 )) + bv(t m + o(t m+1 )) + cw(t n + o(t n+1 )) 
c'(t) = av(t + o(t 2 )) + b'v(t m ' + o(t m ' +1 )) + c'w(t n ' + o(t n ' +1 )) 
c"(t) = a"u(t + o(t 2 )) + b"v(t + o(t 2 )) + c"w(t n " + o(t n " +1 )). 

(3.14) 

This time we estimate the product of the characteristic functions for instance by 

#<5iAi (t)0y 2 A' 2 (t)9s»A% (* ) 

while the determinant can be estimated as above and k is given by k = min(m, m', n, n', n")- 
1 > 1 so that we have actually the same situation as in case lb) upon synchronizing this 
time Ai = A' 2 = A'^ =: A -> 0. 
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As for case II) we observe that the nonvanishing of the functional determinant at 
p implies that the map x ee / e » is actually invertible in a neighbourhood of p by the in- 
verse function theorem. In other words, there is only one point (to,^,^') such that 
%ee'e"(to, t'o, to) = p. Moreover, since the determinant is non-vanishing at p, all three 
edges must be distinct form each other. It follows now from our choice of edges that p 
must be a vertex v = e D e' D e" of 7 in order that he result is non-vanishing and thus 
from the choice of parametrization t = t' Q = tg = 0. 

Therefore, if we take the limit A" — > first in any order then the condition XA"{p, x e e'e") = 
1 will actually imply xa(p, e) = XA'ip, (e + e')/2) = 1 for small enough A'/ so that we 
can take these characteristic functions out of the integral and replace them by 1 if p is 
a common vertex of all three edges. Also we can replace the operator O ee ' e »(t,t' ,t") by 
O ee ' e "(v). This holds only if the triple intersects in p. 

If not all of e, e', e" intersect in p then the limit will vanish anyway if we take a suitable 
limit of the Aj as we have shown before. We can account for that case by replacing 
Xa{p, e), Xa'(p, (e + e')/2) by Xa(p, v )xa'(p, v). Here v is the common vertex at which 
the distinct e, e', e" must be incident otherwise they could not even pass through a small 
enough neighbourhood of p. We can also assume that all three edges have linearly inde- 
pendent tangents at v and expand still around t = 0. The remaining integral divided by 
Ai'Alj'Ag then tends to 

d 3 tdet(^^)5 {3) {p,x ee , e ,,) = s(e,e',e") [ d 3 x5 {3 \p, x) = ~s(e, e', e") (3.15) 



/[o,i] 3 <9t 
where 

s(e, e', e'% := sgn(det(e(0), e'(0), e"(0))). (3.16) 

The factor 1/8 is due to the fact that in the limit A" — ► we obtain an integral over 
C(e, e', e"), the cone based at p and spanned by e(0), e'(0), e"(0) where the orientation is 
taken to be positive. This integral just equals J R 3 d 3 t5(0, t)5(0, t')5(0, t") = 1/8 as one 
can easily check. This furnishes the proof. 
□ 

We conclude that ( |3.6| ) reduces to (in particular, the operators 012,30^2302,310123 drop 
out) 

Jm E(p, A, A', A")/ = g , 3 rJ 1(A) 'J ( ;^ A(P. t ') y O„,,,(0, 0, 0) 

where t> on the right hand side is the intersection point of the triple of edges and it is 
understood that we only sum over such triples of edges which are incident at a common 
vertex. Moreover, 

6 e , eV ,(0,0,0) = ^e ijk Xi„Xi,X k e andX: := X\h e (0, 1) := trfo/^O, 1) — 9 - (3.17) 

is a right invariant vector fieldQ in the Tj direction of su(2), that is, X(hg) = X(h). The 
factor of 1/8 in ( 3.17 ) comes from the fact that at t = t' — t" = we get only one-sided 



functional derivatives of the various edges involved (m(0, 0,0) = 1/8, see above). We 

2 This important observation, which is fundamental for all that follows, is due to Abhay Ashtekar and 
Jurek Lewandowski 
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also have extended the values of the sign function to include which takes care of the 
possibility that one has triples of edges with linearly dependent tangents. 
The final step consists in choosing A = A' and taking the square root of the modulus. 
We replace the sum over all triples incident a common vertex Y^ e e' e" by a sum over all 
vertices followed by a sum over all triples incident at the same vertex J2 v &v(-j) Eene'ne"=ir 
Now, for small enough A and given p, at most one vertex contributes, that is, at most one 
of Xa{v,p) 7^ because all vertices have finite separation. Then we can take the relevant 
Xa(p, v) = Xa(p, v ) 2 out of the square root and take the limit which results in 



V(R), = I d 3 p^det(q)(p) = f d 3 pV(p), 

J R < J R 



veV(-y) 



v v , 



i 3 

£ S ( e ' e '' e ")Qee>e" 



e ijkX l e X J e ,X e „ (3.18) 



d ' ' ° e,e',e"e£;(7),ene'ne"=u 
H v j v k 



where we could switch the order of the X's because a triple contributes only if the corre- 
sponding edges are distinct and so the X's commute. 

Expression ( p.!8| ) is the final expression for the volume operator and coincides up to a 
factor ^27/8 with the expression found in |5|, [| while it is genuinely different from the 
one found in @] as pointed out in [f|. Note that the final expression is manifestly diffeo- 
morphism covariant. Although the procedure of adapting the limiting to a given triple 
of edges is somewhat non-standard there is an argument in favour of such a procedure 
: the discussion in lemma ( |3.1| ) reveals that any other regularization which would result 
in a finite contribution for the case where s(e, e'e") is zero would necessarily depend on 
the higher order intersection characteristics of a triple of edges. However, since such a 
quantity is not diffeomorphism covariant which is inacceptable, the dependence must be 
trivial, that is, a constant (this is precisely what happens in ( |3.4| ) according to fl3.11| ). 



While there is no a priori kinematical reason to prefer one operator over the other, there 
is a dynamical reason : the volume operator can be seen as an essential ingredient in 
the regularization of the Lorentzian Wheeler-DeWitt constraint operator P, [|. The fact 
that ( p. 18 ) contains s(e, e', e") is the precise reason for this operator to be anomaly-free ! 



Therefore, the regularization leading to ( |3.18| ) is the one singled out by the dynamics of 
the theory ! 



4 Spectral Analysis 

In this section we summarize a number of results due to Abhay Ashtekar and Jurek 



Lewandowski obtained in the series of papers [|T^, [|, [|. We have included them here for 
the sake of self-containedness of the present paper. More specifically, the cylindrical con- 
sistency of the family of volume operators V 1 as defined in [[L8[] and in the present paper, 
the self-adjointness of the corresponding projective limit and the fact that the spectrum 
of the volume operator is discrete was first observed in |§ . The proof of these properties 
given in |J uses methods developed in jIS[ 0, §, [|. 



10 



This section is subdivided into three parts. First we prove that the family of opera- 
tors derived in (|3.18 ) defines a linear unbounded operator on 7i. Next we show that the 



operator is symmetric, positive semi-definite and admits self-adjoint extensions and finally 
we show that its spectrum is discrete and that the operator so defined is anomaly-free. 

4.1 Cylindrical Consistency 

What we have obtained in ( [3.18] ) is a family of operators (V(R) 1 , D 7 ) 7e r- That is not 
enough to show that this family of cylindrical projections "comes from" a linear operator 



on H. As proved in fll8[ , for this to be the case we need to check that whenever 7 C 7' 
then 

1) p*yD^ C Dy where p yy > is the restriction from 7' to 7. This condition makes sure 
that the operator defined on bigger graphs can be applied to functions defined on smaller 
graphs. 

2) (V(R)y)\~/ = V(R)~f, this is the condition of cylindrical consistency and says that the 
operator on bigger graphs equals the operator on smaller graphs when restricted to func- 
tions thereon. 

A graph 7 C 7' can be obtained from a bigger graph 7' by a finite series of steps consisting 
of the following basic ones : 

i) remove an edge from 7', 

ii) join two edges e', e", such that e'fle" is a point of analyticity, to a new edge e = e'o(e") _1 

iii) reverse the orientation of an edge. 

Clearly, a dense domain for V(R) 1 is given by D 7 := Cy\^(A/G)- This choice trivially 
satisfies requirement 1) since functions which just do not depend on some arguments or 
only on special combinations h e = h e ih e ",h e i = h~ x are still thrice continuously differen- 
tiable if the original function was (here we have used the fact that 577(2) is a Lie group, 
that is, group multiplication and taking inverses is an analytic map). 
Next, let us check cylindrical consistency. Consider first the case i) that 7 does not depend 
on an edge e on which 7' does. Then clearly 7Q/ 7 = for any function cylindrical with 
respect to 7 and so in the sum over triples over vertices in (|3.18|) the terms involving e 
drop out. 

Next consider the case ii). If e = e' o (e") -1 is an edge of 7 and e', e" are edges of 7' where 
v := e' PI e" is a point of analyticity for 7 while for 7' it is not, then while v is a vertex 
for 7' it is only a pseudo-vertex for 7 and so in V{D) 1 there is no term corresponding to 
v. On the other hand, since the vertex v is a pseudo vertex for 7 it is in particular only 
two-valent and so the corresponding term in V(D)y drops out. Likewise, if v is a vertex 
for 7 at which the outgoing edge e is incident, then from right invar iance of the vector 
field we have X e = X e / f e n\-i = X e i and so at vertices that belong to both 7 and 7' the 
corresponding vertex operators coincide. 

Finally, case iii) is actually excluded by our unambiguous choice of orientation. 

We conclude that there exists an operator (V(R),D) on TC which is densely defined on 

D = Cj\\A/g). 

4.2 Symmetry, Positivity and Self-Adjointness 



Notice that the vector field iX e is symmetric on 7i 7 , the completion of Cyl (A/Q) with 
respect to ^0,7; e an edge of 7, because the Haar measure is right invariant. It follows 
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from the explicit expression ( |3.18| ) in terms of the iX e that all the projections V(R)^ 
are symmetric. In this special case (namely, the volume operator leaves the space _D 7 
invariant) this is enough to show that V(R) is symmetric on D. 

Furthermore, all V(i?) 7 are positive semi-definite by inspection so that V(R),D is a 
densely defined, positive semidefinite and symmetric operator. It follows that it has self- 
adjoint extensions, for instance its Friedrich extension. 



4.3 Discreteness and Anomaly-freeness 



The operator V(D) has the important property that it leaves the dense subset Cyl^ ) (A./Q) C 
Ji invariant, separately for each 7 G T. Moreover, recall that TC possesses a convenient 
basis, called the spin-network basis in the sequel |21], |22], |23 . 



Let us recall the most important features of that basis. Given a graph 7, let us associate 
with each edge e G E(^y) an spin quantum number j e > 0. Also, with each vertex v G ^(7) 
we associate a certain contraction matrix c v as follows : if ttj is the irreducible represen- 
tation of SU(2) corresponding to spin j then form the tensor product (S>eG£(7) 7r j e (^e) and 
for each vertex v contract the indices A of the matrices 7Tj e (h e )AB for all e incident at 
v in a gauge invariant fashion. We obtain a gauge invariant function T^j s (A), called 

a spin-network function. It turns out that given 7,j there are only a finite number of 
linearly independent c compatible with 7,j. Moreover, one can show that spin- network 
functions defined on different graphs are orthogonal and even on the same graph they 
are orthogonal whenever the vectors j don't equal each other. It follows that there is an 
orthonormal basis 

< T 7J > T yJv >= S jn' 5 U'ks' ■ 

Now it is obvious that the operator V(R) leaves the finite dimensional vector space U 
spanned by spin-network states compatible with j,j invariant. The matrix 

(v(R),j)c,;> ■=< t,j, s \v(R)Ki, > (4.1) 

is therefore finite-dimensional, positive semi-definite and symmetric. The task of comput- 
ing its eigenvalues therefore becomes a problem in linear algebra ! 
Next, since from ( |3 .18 ) 

V(R), = £ 3 P £ K l7 (4.2) 

t>eV(7)ni? 

and since V vn involves only those e G E{pf) with v G e, we find that V vri can only change 
the entry c v in c. In other words, [V v ,-y, = and each V Vjl can be diagonalized 
separately. 

Finally, since the spins j e only take discrete values it follows that 7i 7 has a countable basis 
and the spectrum that V(R) attains on _D 7 is therefore pure point. Let us check whether 
this is the complete spectrum. Assume it were not and let P the spectral projection on 
the rest of the spectrum (the existence of the spectral projections relies on the fact that 
V(R) is self-adjoint and not only symmetric). It follows that u = Pv is orthogonal to 
where v is any vector in 7i 1 . But _D 7 is dense in 7i 7 and so we find for every e > a <p G D 1 
with I \u — 0|| < e. Now we have from orthogonality e 2 > | \u — <p\ \ 2 = \ \u\ \ 2 + \ \<p\ | 2 > | \u\ | 2 
and so u = 0. This shows that the complete spectrum is already attained on D 1 . It is 
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purely discrete as wellQ. 

Last, we wish to show that the volume operators are anomaly free (given the fact that we 
have largely adapted our regularization to a graph, this statement is far from trivial). By 
this we mean the following : given any two open sets R\, R2 C £ we have vanishing Poisson 
brackets {V (R\) ,V (R2)} = because the functionals V(R) depend on the momentum 
variable E^(x) only. Now, given a function / cylindrical with respect to a graph 7, it 
is not at all obvious any more that [V(Ri), V(R 2 )]f = for any such /. Fortunately, 
given the above characterization of the spectrum, the commutator can be easily proved 
to vanish on cylindrical functions. To see this, note that the above results imply that 
if we choose any region RIj) such that 7 C R(j) then there exists an eigenbasis of 7^ 7 
of V{R{ r f)). Now consider any region R. Since all regions are open by construction, all 
regions fall into equivalence classes with respect to 7 : R, R! are equivalent if they contain 
the same vertices of 7 (any vertex either has a neighbourhood which lies completely inside 
R or it lies outside). Therefore any two V(R),V(R') differ at most by some of the V v ^ 
all of which are contained in the expression for V(R( / y)). Since the V vn commute the 
eigenbasis of V(R{pj) is a simultaneous eigenbasis of all V va for all v G ^(7) and so this 
eigenbasis is a simultaneous eigenbasis of all V(R) T Since all 7i 7 are orthogonal, we have 
a simultaneous eigenbasis for all V(R). 

While it is in general not enough to verify that two self-adjoint, unbounded operators 
commute on a dense domain (rather, by definition, we have to check that the associated 
spectral projections commute) in our case we are done because the spectral projections 
are the projections on the various D 1 because the point spectrum is already the complete 
spectrum. Thus we have verified that the commutator algebra mirrors the classical Poisson 
algebra. 



5 The complete set of matrix elements 



In this final section we are going to compute ( |4.1| ) in a spin-network basis for each 
7, j, c, c!, R. The advantage of our approach as compared to the ones proposed already 
in the literature JTOf will be obvious : in one works with an overcomplete and non- 



orthogonal set of states. This prevents one from giving a unified treatment as it is done 
here and only allowed one to find a formula for the eigenvalues in a few simple cases. 
More generally one has to check all the time the linear independence of the states and 
can compute the matrix elements only in a case by case analysis. 

In the next subsection we show in more detail a result which follows immediately from 
0, [5], |H : namely, that the problem of computing the matrix elements of the volume 
operator reduces to that of computing the matrix elements of a homogenous polynomial 
of degree three of spin operators for a spin system. 

This opens access to powerful techniques well-known from the quantum theory of 
angular momentum for spin systems. Then we will show how the actual computation 



3 in the physical sense that it is attained on a countable basis so that the eigenvalues only comprise a 
countable set. In a mathematical sense one would need to check that there are no accumulation points 
and no eigenvalues of infinite multiplicity. To settle this question it would be enough to show that the 
formula for the eigenvalue is an unbounded function of j :— J^eesfa) 3 e because every countably infinite 
set of j's corresponding to different choices of j diverges. This is one possible future application of the 
explicit matrix element formulae which we derive in the next section. 
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can be done by viewing the problem as a task in spin recoupling theory and that all the 
matrix elements can be written as algebraic functions of 6j symbols for which a closed 
formula exists (Racah formula). Finally we will diagonalize the volume operator for an 
infinite number of (very special) graphs. 

5.1 The left regular representation and the spin- network repre- 
sentation 

Consider a spin- network state T 7 jg(A). This is a state for the connection representation 

L2(-A/G, dfio) in the sense that the (equivalence classes modulo gauge transformations of) 
connection operators are diagonal in this representation. Notice that a spin-network state 
is gauge invariant so that it can be viewed alternatively as a function of [A], the gauge 
equivalence class of A. Here we will use the fact that df^o extends to A pOJ and reduces 
to the measure indicated in the introduction when integrating gauge invariant functions. 
Consider the usual Dirac (generalized) states 5a defined by Sa(A') := S(i (A, A) where <5 W 
is the 5 distribution with respect to /i P3[ and A, A e A. We define the spin-network 



representation by 

<U(T tJi? ) := [ d^ (A')5^(A)T lls (A') = T^(A) =:< A| 7 , j, c > . (5.1) 

So far this is only notation. The advantage of this notion will become transparent upon 
proving that the states \ c> are nothing else than the usual angular momentum states 
associated with an abstract spin system. 

To see this, let v be a vertex of 7 and ei,..,e n the edges of 7 incident at v and we 
have chosen orientations such that all of them are outgoing at v. The spin-network state 
T-,?JA) can be written 

n 

T~fj } g(A) = (C 1 ,) mi .. mn J^(7T Je . (/l e4 (yl))) m - m ^(M 1 ,) m ' i .. m / i 

8=1 

= tr(c v ■ ®^ 1 7r jei (h ei (A)) ■ M v ) (5.2) 

where c v is the vertex contractor corresponding to v which contracts the indices corre- 
sponding to the starting points of the ej and M v contracts the indices corresponding to 
the endpoints (recall that the holonomy is a path ordered exponential with the smallest 
parameter value to the left). Note that c v is C— valued while M v still depends on the rest 
of the graph. 

Consider a local gauge transformation at v, that is, h e . — > gh ei for some g G SU(2). Then 
( |5.2|) becomes ti(c v ■ <S>" =1 7rj e . (g) ■ ®i = i^j e .(h ei (A)) ■ M v ). In order for this to be gauge in- 
variant, notice that the tensor product representation ®iTij e _ (g) can be written as a direct 
sum of orthogonal irreducible representations iij(g), and so we just need to choose c v to 
be proportional to one of the various equivalent (but orthogonal) trivial representations 
^0(9) = ^(l) contained in that decomposition (see |23| for more details) in order for (|5.2|) 
to be gauge invariant. This is precisely what a spin-network state is. 
The point is now that 

n 

Il( 7r ie i (^ei(^)))m i m< =:< A\(rrii, ji), .., (m n , j n ); mi, .., m' n > (5.3) 
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transforms precisely as the usual angular momentum eigenstates 



|(mi,ji),.., (m^jVOjfci., ..,k m >:= Y[\mi,ji;ki,..,k m > 



i=l 



of an abstract spin system of n spins ji,..,j n under rotations (see, e.g., pMfl). This is 
because \(m,j); k > transforms as U(g)\(m, j); k >= iij(g) m/m >\(m' , j); k > by definition, 
where U(g) is a unitary representation of the rotation group. Here ki,..,k m are some 
additional quantum numbers of observables commuting with the operators Jj which in 
our case coincide with the m^, ..,m' n . In other words, the left regular representation 
L g f(h e ) = f(gh e ) that we are dealing with in the connection representation can be seen 
to be equivalent with the representation U(g) of the spin-network representation. 
Let us now couple the spins ji, ..,j n to resulting spin j = and denote a particular state 
with zero angular momentum by |0 >. We claim that up to a constant 

=< (mi, ji), .., (m n ,j n ); fa, .., k m \0 > (5.4) 

where the inner product is now to be understood in the angular momentum Hilbert space 
of the abstract spin system. In other words, the c„ are just usual Clebsh-Gordan coeffi- 
cients. The proof is easy, we have 

{Cv)m 1 ..rn n IIj=l ^j; (sOrr^m^ (<?) 

= < ^(flO(( m 'l, jl), ■; On, 3n)] fa, .., k m )\0 > 

=< (m' 1 ,ji),..,(m' ri ,j n );/ci,..,A; m |f>(^- 1 )0 > 

{('-v)m',..m' n 

where we have used unitarity of the group and of the representation and that U(g)\0 >= 
|0 > is rotation invariant. The factor of proportionality is seen to equal unity because 
c 2 v = 1 and since E mi .. mn \{n r ii,ji)..{m n ,j n );k 1 ..k m >< {m 1 ,j 1 )..(m n ,j n );k 1 ..k m \ = 1 due 
to completeness and Wigner-Eckart theorem [|J] . 

Next, consider the self-adjoint right invariant vector fields Y e := —i/2X e . We have 
[X 1 , X- 7 '] = —2e^ k X k and so = %^ k Y k which is the usual angular momentum 

algebra. It follows that (we use a dual interpretation of the SU(2) indices A,B,C, .. 
which we may choose to take values ±1/2 and so can be interpreted to be the eigenvalues 
of the 3-component of angular momentum) 

\j = 1/2, C > Yjh e (A)cD =< 3 = l/2,C\(-i/2)(rMA))cD 
= (i/2)<( Tj ) EC (l/2,C)\h e (A)h D 

= (i/2)j t . „<U(eMtTj)(l/2,E)\h e (A)) ED =< \j2,E\J 3 h e {A) ED (5.5) 



t=o 



and so by multilinearity of the tensor product representation we find due to self-adjointness 
of Ji (set Yj :=Y ei ) 

id - -> 

= < j v ,m >< A\j,j v ,m; m\JjO > (M v ) m , v .. m , n (5.6) 

where j v comprises the components of j which correspond to the edges incident at v. 
Since spin-network states are thus characterized by the various zero-angular momentum 
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eigenstates |0 > of the abstract spin system, there is a one-to-one correspondence between 
the possible choices of c v and |0 >. So let |0 > correspond to c v and |0' > to c' v . Let p(Y) 
be any polynomial of the Yj corresponding to the vertex v then we have for its matrix 
elements 

<T J>ls \p(Y)\T^> 

= / d[Mo{A)<j v ,m\Q >< A\jJ v ,rh;rh' > (M v ) m > m , <j v ,m\p(J)0' > x 

JA/g 1 

x < A\j,j v ,m; rh > (M v ) 

m' 1 ,..m' rl 

= < j v , rh\0 > < j v , m\p(J)ty > n 5 c v ,,c' v , 

= < ob(j)t ' > n K'K' =< p^°\°' > n k><, ■ ( 5 - 7 ) 

Here we have used ( |5.6|) in the first step, the orthonormality of the spin- network states in 
the second step and in the last step we have again recalled that J2m \jv,fh >< jv, ™>\ = 1 
is the identity on the abstract angular momentum Hilbert space labelled by the spins j v 
(completeness, notice that we sum over rh in (|5.7p ). 

Expression ( |5.7| ) is a huge simplification : in order to compute the matrix elements of the 
volume operator it is enough to compute the matrix elements of the angular momentum 
operator polynomial eijkJ % e J 3 e iJ^n between states that have zero angular momentum coming 
from j v . This task is devoted to the next subsection. 



5.2 Recoupling Theory for n angular momenta 

The first task is to label all the linearly independent states |0 > which correspond to 
the coupling of n angular momenta ji, ..,j n to zero angular momentum. This is a well 
known problem in the quantum theory of angular momentum and runs under the name 
recoupling theory. We review here what we need for our computation. 
The states \ j, rh > are a complete basis for a spin system consisting of n angular momenta. 
In this basis the 2n mutually commuting operators (J)) 2 , Jf are diagonal with eigenvalues 
+ l),mj. 

We are obviously interested in a basis in which the square of the operator of total angular 
momentum J := J\ + .. + J n and its 3— component are diagonal with eigenvalues j(j + 
l),m. In order to label that basis we need 2(n — 1) more quantum numbers of mutually 
commuting operators and commuting with J 2 , J 3 . The n operators Jf are readily verified 
to satisfy these conditions but we need n — 2 more. In order to motivate our choice notice 
that V vn is the square root of the modulus of i(3/2) 3 ^ J2i<j<k e ( e i, e J, ^k)qijk where 

-^Iijk '■= gtijkX^X^Xx = [YjY},YjY^] = -[(F/j) 2 , (Yj K ) 2 ] 

where Yjj := Yi + Yj. Here we have used total antisymmetry of e(ei,ej,eK),qijK and 
the fact that [X/, Xj\ = 0, 1 ^ J in order to get rid of 3! and in order to sum only over 
I < J < K. This observation^ urges us to study the operator 

qijK ■= [Jui j jk] where J u = J/ + Jj (5.8) 

4 This fact that the third order polynomial tv{X j[X j 1 X k]) in the X can be written as a commutator 
of two second order polynomials tv(XiXj) which are gauge invariant was made independently also by 
Abhay Ashtckar and Jurek Lewandowski 
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and thus it would be convenient if we could work in a basis in which the Jjj also were 
diagonal. Clearly we cannot achieve this for all n(n— l)/2 operators Ju but we can work 
in different bases adapted to <Juk for a particular triple I < J < K and expand them in 
terms of each other. 

We will develop a theory of recoupling schemes here which is somewhat different from 
the one encountered in the standard literature and we also consider arbitrary n. The 
following paragraphs will therefore be somewhat detailed. 

Definition 5.1 A coupling scheme based on a pair 1 < I < J < n is an orthonormal 
basis \g(IJ),j,j,m >, diagonalizing besides J 2 , J3, Jj, I = l,..,n the squares of the 
additional n — 2 operators G 2 , ■ ■, G n _\ where : 

G\ :— Ji, G2 '■= G\ + Jj, G3 := G2 + Ji, G4 := G3 + J2, ■ ■, Gj+i := Gj + J/_i, Gj +2 '■ = 
Gi + i + J1+1, Gj +3 := Gi +2 + J1+2, ■■, Gj := Gj-i + Jj-i, Gj + i := Gj + Jj+i, Gj + 2 '■— 
Gj+i + Jj+2, ■ ■, G n := G n -\ + J n = J. 

The vector g(U) :— (g2(U), ■■, g n -i{U)) denotes the eigenvalues of the squares of G2, ..,G„_i. 
We will call the coupling scheme based on the pair (12) the standard basis. 

Let us check that the Gj satisfy the angular momentum algebra and that the squares of 
the operators J, Ji, Gi and the operator J 3 are mutually commuting : for this it will be 
enough to check that the operators for the standard basis are mutually commuting since 
there is a relabelling of indices such that any coupling scheme becomes the standard basis. 
Since [Jj, Jj] = 0, / ^ J we find [G\, G\] = Ej=i[4, J j) = i<Hjk Ej=i J J = ^kG*} which 
implies that G\ has spectrum gi(gi + 1) with gj integral or half integral. 
Let I < J, then we can write Gj = Gj + + •• + Jj) —: Gj + G\ so that [Gj, G'j\ = 
because G^G'j involve disjoint sets of J K 's and so [Gj,Gf\ = [Gj + [G'jY^G 2 ] + 
2(G' I ) i {G j I [G I , G{] + [G), G\}G{} = 2i{G' I ) i t ljk {G\G k I + G^Gj} = 0. The same argu- 
ment can be used for the other commutators. 

The virtue of this definition is the following : we have now a neat labelling of the vec- 
tors |0 > referred to in the previous subsection : we can choose them to be the vec- 
tors of the standard basis corresponding to j = 0, that is, they are just the vectors 
\g(12),j,j = 0,m = > (m = is forced by j = 0) and the different vertex contractors 
just correspond to the various possible values of the intermediate (or so-called re-coupling) 
vectors G 2 , .., Gn-i- We will therefore compute the matrix elements of qjjK in this basis. 

— * 

In the sequel we will drop the indices j,j = 0,m = since they are the same in every 
coupling scheme and since qjjK commutes with J 2 , Jj,(J) 3 matrix elements between vec- 
tors with different values of these operators therefore vanish anyway. We will also write 
g instead of g(l2). 

After these long preparations we are now in the position to complete the derivation. We 
use completeness of the bases \g(IJ),j,j,m > and have 

< 9\<iiJK\f? > 

= ]T g 2 (IJ)(g2(IJ) + 1)[< g\g{U) >< (9ij)\JjkW >~< MkWu >< 9(U)\g >] 
]T g 2 (IJ)(g 2 (IJ) + l)g 2 (JK)(g 2 (JK) + 1) < g{IJ)\g" >< g{JK)\g" > x 

g{IJ),g(JK),g" 

x[< g(U)\g>< g{JK)\$ >-< g(JK)\g>< g{IJ)\$ >] (5.9) 

where we have used that the so-called 3nj— symbols < g(IJ)\g(JK) > are real (also 
< g(IJ)\g(KL) > for any I < J, K < L are real). To see this, insert a unit in terms of 
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the basis \j,rh > which allows us to write the 3nj— symbol in terms of Clebsh-Gordan 
coefficients which are real [24}]. It is understood that in (|5.9| ) and in what follows the 



summation extends only over the allowed values of g etc. which are purely fixed by the 
values of ji, ..,j n . 

The task left is to write the 3nj— symbols of the form < g(IJ)\g' > in terms of known 
quantities. We will derive below a formula which expresses it in terms of the 6j symbol 
for which a closed formula (the Racah formula P^j ) exists. 
We will denote g(IJ) = (g 2 (ji,jj), gs(g 2 , ji), ••, 9i+i(9i,ji-i), 91+2(91+1, ji+i), 
gj(gj-i,3j-i),9j+i(gj,jj+i), »g n -i{g n -2, 3n-i)) where the notation means that we 

couple J, J' to resulting spin j" of J" thereby keeping explicit track of the coupling scheme. 
The following two elementary lemmata will be crucial in computing the 3nj— symbols. 

Lemma 5.1 < g(IJ)\g' > 

=< 92 (ji, j ; j) , 93 (92, ji) , ■ ■ , gi+\ (gi, ji-i) , gi+2 (gi+i , ji+i ) , ■ ■ , g j (gj-i,jj-i)\ 
\g' 2 (Juj2),g 3 (g 2 ,j3), -,gj(g'j-iJj) > s gj,g' J - S 9n-i,g' n _ 1 where on the ri 9 ht hand Slde we con - 

sider a spin system consisting of spins J±, .., Jj and total spin quantum number gj. 
Proof : 

Since Gj(12) 2 , .., G n _x(12) 2 are diagonal on both vectors it is immediately clear that the 
3nj— symbol must be proportional to the Kronecker delta 8 gKtg i for K = J, ..,n — 1. To 

get the factor of proportionality we expand both vectors in terms of the basis \ j, m > and 
Clebsh-Gordan coefficients. To do this, note that by definition 

< 02 (ji, 32), ■-, g n -2{g n -3, jn-2); g n -im - m n ,j n m n 

\g2UlJ2), ..,g n -2(gn-3,jn-2),gn-l(gn-2,jn-l),j(gn-l,jn),m > 

=< gn-i(g n -2, j n -i)m — m n, j n m n \j, m; g n -i, j n > where on the right hand side we have 
the standard Clebsh-Gordan coefficient. Here 

\g2U1J2), -,gn-2{g n— 3; jn— 2)) 5Vi— 1^ TTlni jvX^n ^ 

:= Ifl^O'i) h), ■■,gn-2(gn~3,jn~2),gn-i(gn-2,jn-i),m- m„ > ®|j„m n >. We can now iter- 
ate this procedure and obtain 

WzUhh): ■■> g'n-2(g'n-3, jn-2) , g n -l(g'n-2, jn-l) , j (g'n-!, jn) , TR > 

= Y,mj+..+m n = m < g' n -i m - m n,jnm n \jm >< g' n _ 2 m-m n -m n _i , jn-xm^-i \g' n _ x m-m n > 
•• <g'j'mj,jj +1 mj +1 \g'j +1 ,Tnj + mj +1 > 

\g2Uuj2), ■■,gj(g'j-i,jj),mj > ®\j J+1 ,m J+ i > ®.. <g> \ j n ,m n >. Similarily we have 

\g2(jl,jj), ..,g n -2(gn-3,jn-2),gn-l(gn-2jn-l),j(gn-l,jn),m > 
= Tlmj+..+rn n =m < 9n—lTR> ~ m n, jnTnn\j m > 

< 5 , n-2^-^n-^n-l,jn-l"ln-l|5'n-l^-^n > •• <gjmj,jj +1 mj +1 \gj +1 ,mj + m J+1 > 

\g2{ji,jj),-;gj{gj-i,jj-i),mj > ®\j J+1 ,m J+1 > ®..®\j n ,m n >. Thus we have the inner 
product 

< g(JJ)\g > = J2mj+..+m n =m 

< g2(jijj),--,gj(gj-ijj-i),rnj\g 2 (ji,j 2 ),--,gj(gj-ijj),rnj > x 

x < g' n _ x m - m n ,j n m n \jm >< g' n _ 2 m - m n - m n _i, jn-\m n -i\g' n _ x m - m n > .. < 
g'jmjjj +1 mj +1 \g' j+v mj + mj +1 > x 

x < g n _ x m - m n ,j n m n \jm >< g n - 2 m - m n - m n _i, j n _im n _i|p n _im - m n > .. < 
9jmj,jj +1 m J+1 \gj + i,mj + m J+1 > 

=< g2(jijj),--,gj(gj-ijj-i),m\g 2 (jij 2 ),.. } g'j(g , j _ 1 jj),m> x 

x E mj+1 ,..,m„ < g' n -i m ~ m n , j n m n \jm >< g' n _ 2 m -m n - m n _ u jn-irrin-^g'^m - m n > 
•• <g'jmjJj +1 mj +1 \g' J+1 ,mj + m J+1 > x 
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x < g n _im - m n ,j n m n \jm >< g n - 2 m - m n - m n _i, j n -im n -i\g n -im - m n > .. < 
gjmj,jj +1 mj +1 \g J+1 ,mj + m J+1 > 

=< g2(jiJj),-;gj(gj-ujj-i),rh\g' 2 (j u j2),.. J g'j(g' J _ 1 Jj),m > 8 gjiglj ..5g n _ ugln l where in 
the second step we have used the fact that the 3nj— symbols are independent of m 
(Wigner-Eckart theorem^]) so that m is arbitrary and in the last step we have used the sum 
rule for the Clebsh-Gordan coefficients, that is, J2m 1 +m 2 =m( < 3i m u jif^i^j^i >) 2 = 1. 
□ 

Lemma 5.2 < g 2 (j u j 2 ), .., gxiSK-u 3k), 9k+i(9k, 3k+i, 9k+2(9k+i, 3K+2);m\ 

\92(jlj2),-,9K{9K-l,jK), 9k+1 (9K,jK+2),g' K+ 2(9'K+ldK+l)^ m > 

=< gK+i(gKjK+i),gK+2(gK+ijK+2);m\g' K+1 (g K ,j K+2 ),gK+2(g'K+iJ3)i m > 

Proof : 

Again we expand 

\92, ■■,gK+2] m>= J2m K+1 ,m K+2 

< g K m-m K+l -m K+ 2,3K+im K+ i\g K+ im-m K+ 2 >< gK+2m-m K +2, 3K+2m K +2\gK+2m- 
m K +2 > x 

x\g 2 , ..,g k ;m- m K+1 - m K+2 > ®\jK+im K +i, 3K+2m K+2 > and 

\g2,-,gK, 9k+V 9K+2', m>= J2m K+1 ,m K+2 

< 9Km-m K+1 -m K+ 2,jK+2rn K+ 2\g'K+i m - m K+i >< 9K+2m-m K+1 , j K+1 m K+1 \g K+2 m- 
m K +i > x 

x | <7 2 , ••, 9k\ m ~ mK+i — fnR+2 > <%>\jK+i'm>K+i, 3k+iTtik+2 >■ We find for the inner product 
of these states due to the normalization just 

J2m K+1 ,m K+2 < 9Km - m K+1 - m K+ 2,j K +im K+1 \g K+1 m - m K+2 > x 
x < g K+2 m - m K +2, jK+2m K+2 \g K+2 m - m K+2 > x 

x < g K m-m K+ i-m K+ 2 } jK+2rn K+ 2\g' K+1 m-m K+ i >< g K+2 m-m K+1 , j K+1 m K+ i\g K+2 m- 
m K +i > 

=< 9k+i(9k, jK+i), 9k+2(9k+i, jK+2); m\g' K+1 (g K , j K+2 ), g K+2 (g' K+1 , j 3 ); m > as claimed. 
□ 

We are now ready to reduce out the general 3nj— symbol. We assume that gx = g'x, K = 
J, ..,n — l as otherwise we get zero by lemma 5.1 . Upon repeatedly using the complete- 
ness relations of the coupling scheme bases we first shift jj and then jj to the right in the 
matrix element 

<g(U)J> 

=5.1 < gb(ji, 3 j), 93(92, ji), ..,gi + i{gi,ji-i),gi+ 2 {gi+i,ji+i), ~,gj{gj-i,jj-i)\ 

I g'2 tii , h ) , 9 3 (92 » h ) > • • , 9 j- 1 (9 'j-2 > 3 ] j- 1 ) , 9 j (9 j- 1 , 3 \j ) > 
=0 12 <92UiJj)>9a(92ji),g4,--,9j\h2{3i>h)>93(fo,jj),9t,",9J > 

h 2 

< h 2 {ji,ji),gz{h 2 ,jj),g^ .., gj\g' 2 {ji, 32), gM^h), -,g'j-i,gj > 
=0 12 < 92(31, jj), 93(92, ji)\h 2 (ji, 3i),93(h 2 , jj) > x 

h 2 



5 Proof : Let us expand \g(U) >— J2g(KL) < 9{KL)] j,m\g(I J)] j,m > \g{KL)\j, m > where we 
assume that m < j w.l.g. Now apply the ladder operator J' to both sides of this equation. We find 
J' \g{MN) ;j,m >= c(j, m)\g(M N) ; j, m + 1 > for any K < L where the coefficient only depends on j, m. 
Comparing coefficients we find < g(KL);j, m\g(IJ);j, m >—< g(KL); j,m + l\g(IJ);j, m + 1 >, that is, 
the 3nj— symbols are rotation invariant. □ 
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5.1 



x < hzUiJi), 9s(h 2 , jj), 94(93, 32), 9s, •■,9j\h2tii,ji),h3(h2,j2),94(h 3 ,jj),gB, ~,9j > 

< h 2 (ji, ji),h 3 (h 2 , j2), g4(h 3 Jj), gs, ■■, gj\g 2 (ji, h), gM, h), ■■, g'j-i, 9j > 
§ E< 92(jiJj)>9s(92,ji)\h2(jiJi)>9a(h2jj) > x 

hi 

X E< g3(h2,jj),g4(g3,j2)\h-}(h 2 ,j2),g4(h3,jj) > X 

x < h 2 (jiji), h 3 (h 2 J 2 ), g4(h 3 , jj), g 5 , ■-, gj\g' 2 (ji, 32), g 3 (g' 2 , 33), ■■, g'j-i, gj > 
Z < 92(ji, jj), 93(92, ji)\h 2 (ji, 3i),g3(h 2 , jj) > x 

hi 

x Z < g3(h2,jj),g4(g3,32)\h 3 (h 2 J 2 ),g 4 (h 3} jj) > .. 

h 3 

Z < 9i(hi-i, jj), 91+1(91, ji-i)\hi(hi-i, ji-i),gi+i(hi, jj) > x 

fez 

x £ < 9i+i(hi,jj), 91+2(91+1, 3i+x)\hi+x(h I: 3i+x),gi+2(hi+i,jj) > x 
fei+i 

X Z < gi+2(hl+l, jj), gi+3(gi+2, 3l+2)\hl+2(hl + l, jl+2), gi+3(hl+2, jj) > ■■ 

hi+2 

< 9J-i( h J-2,3j),gj(gj-i,3J-i)\hj-i(hj- 2 ,jj-i),gj(hj-i,jj) > x 

fej-l 

x < h 2 (ji,ji)h(h 2 ,j 2 ), ..,h I (h I - 1 ,j I - 1 ),h I+1 (hi,ji +1 ), ..,hj- X (hj-2,3J-x),gj(hj-i,3j)\ 

\g' 2 (3\, 32), ■■, g'i(g'i-i, 31), g'i+i(g'i, 31+1), ■■, g'j(g'j-i, 3j) > ■ (5.io) 



Now the last matrix element in the last line of ( J5.10 ) is different from zero only for 



Hk = g'j+,K = I,--,J — 1 by lemma |5.1| (we have gj = g'j already) and in that case 
reduces by the same reasoning to 

< h2(jiJi)h 3 (h2j2),--,h I - 1 (h I -2,ji-2),gKh I _ 1 J I _ 1 )\g 2 (j 1 J 2 ) J ..,g' I (g' I _ 1} j I ) > 

= < h 2(jl,jl),h 3 (h2,j2)\k2(jl,j2),h 3 (k2,jl) > x 

kl 

X J2 < h 3(h, ji), hi(h 3 , j 3 )\k 3 (k 2 , j 3 ), h(h, ji) > •• 

J2 < fej-i(*j-2,jj),/ij(^/-i,jj-i)|A;j-i(A:/_2,jj-i),^(A;j_i,ji) > x 
fcj_i 

x < k 2 (3i, 32),h(h, 33), -,^1-1(^-2, 31-1), gi(ki-i, ji)\g 2 (ji, 32), gMJs), ■■, g'i(g'i-i, 31) > 

(5.11) 

and the last matrix element in the last line of ( |5.11| ) is just given by &k K ,g' so that 

the summation drops out. So, altogether, there are only the 1 — 3 summation variables 
h>2, ■■, hi-i. 

Equations ( |5.10|) , (|5.11|) demonstrate that we manage to reduce the 3nj— symbol under 



investigation to a polynomial of matrix elements of the structure 

Type I : < J12 ( Ji , J2 ) , J (j'12 , Js ) I J13 (ji , Js ) , J (j 13 , J2 ) > or 

Type II : < j 12 (ji, 32), 3(312, 33)^23(32, js), 3(323, ji) > • (5.12) 
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In fact, all matrix elements in ( 5 . 1 0| ) are of type I and all matrix elements in (|5 . 11 ) are of 
type I except for the first factor after the first equality which is of type II. 
The point is now that both matrix elements are related to the 6j— symbol |[24|| . The 
6j— symbol is defined through the type II matrix element, namely 

< jl2 (jl , 32 ) , j (jl2 , 33 ) 1 323 (32 ,js), 3 (j23 , 3l ) > 



(2j 12 + l)(2j 23 + l)(-ir^ + ^<j *} J ; f \ . (5.13) 

[ 33 3 323 J 

In order to relate the type I matrix element to the 6j— symbol we need to use the following 
identity for the CG-coefficients 

< jimi,j 2 m 2 \jm >= (-1) J1+J2 ~ J < j 2 m 2 ,jim 1 \jm > (5.14) 

which is a consequence of the choice of phases for the CG coefficients in order to have them 
real valued f24|. Let us now set j[ := j 2 ,f 2 := j u j 3 := j 3} j' 12 := jn,j' 23 := j r3 ,j' ■= j 
then we get 

< ji2(ji, 32), j(ji2, 33)\ji3(ju js), 3(3i3, 32) >=< j[ 2 (j 2 , jtiJ'Uu'faMnUn fatifUnifi) > 

which is almost the standard form except for the wrong order of j[,j' 2 in the first entry. 
Using the expansion 

Em 1+ m 2+ m s=m > < f 2 m 2 , 3[m 1 \j' 12 m 1 +m 2 >< j' 12 m x + m 2 , j! 3 m 3 \fm' > \j[m u .., f 3 m' 3 > 
and ( |5.14j ) we therefore find 

< 3 12 (Ji 2 , j 3 )\j' 23 (j 2, (323, ji) > 

= (-lyi+ii-i™ < j' 12 (j[, f 2 ), f(j[ 2 , f 3 )\j 23 (j 2 , f 3 ), j'(f 23 , j[) > and so we have for the type 
I matrix element 



< 3l2(jl, 32), j(jl2, J3)|jl3 (jl, h), 3(313, 32) > 



(-l) il+i2 - jl2 ^(2ji 2 + l)(2ji3 + 



32 3i 3i2 

33 3 3i3 



(5.15) 



Finally, for the benefit of the reader we note the Racah formula [24 



Jl 32 312 
33 3 323 



A ( ji , 32 , 312) A( 31 , j, j 23 )A (j 3 , j 2 , j 23 ) A (j 3 , j, ,j 12 )w 



3i 32 312 
h j J23 



A(a, b, c] 



\ 



(a + b - c)\(a - b + c)!(-a + b + c)\ 



w 



3i 32 312 



(a + b + c+l)\ 
£(-!)> + l)!x 



, 33 3 323 

X [(n ~ jl ~ J2 ~ jl2)K n ~ jl ~ j ~ j23)K n - J3 - h ~ J23)!( n _ ^3 ~j~ jl2)!] _1 X 

x [(Ji + 32 + J3 + j ~ n)\(j 2 + j'12 + j + j 2 3 - n)\(ji 2 + ji + j 23 + j 3 - n)!] _1 (5.16) 

and the sum extends over all positive integers such that no factorial in the denominator 
has a negative argument. In conclusion, formulae ( |5.9|) , (|5.10|) , (p.llD provide us with 



21 



the general expression for the matrix elements of the operator <Juk in the standard basis 
where it is understood that all occurring matrix elements are known in terms of j\, ..,j n 
via G57L3D , ( |5TT5D , (plSQ . 



Let us call M v jjx(g, g') '■—< glQuxlg 1 > then we have the matrix 

M v (g,g') = ^fj £ e(I,J,K)M vJJK (gJ) (5.17) 

A 4 /<J</sT 



which by inspection of (|5.9|) is a Hermitean matrix of a special kind, namely, it is of 
the form of i times a real-valued skew matrix. It's eigenvalues therefore come in pairs 
±A, A > or are zero. The eigenvalues of the positive semidefmite matrix M^M V are 
therefore |A| 2 or where A is an eigenvalue of M v . We conclude that the eigenvalues of 
K,,7 are given by ^\J\\\ or where both eigenvectors of M v corresponding to eigenvalues 

±A produce the same eigenvalue ^y|A| of V va . We see that the problem of diagonalizing 
the volume operator becomes equivalent to diagonalizing the matrices M v which can be 
done for each v G V(~f) separately. 

In the next subsection we will derive the spectrum in the most simple cases. 
Remark : 

We have nowhere made use of the fact that we are only interested in gauge invariant 
functions. Therefore, our results carry over, word by word, to the case of non-gauge 
invariant functions provided we use the notion of an extended spin-network state. 

Definition 5.2 An extended spin-network state is a state of the form j s as before just 
that the vertex contractors now are Clebsh-Gordan coefficients of the form 
< j 1 m 1 , .., j n m n \j u .., j n , g 2 , ..,g n -i,j > where j maybe different from zero. 

This is important for applications in quantum gravity [Q, [3| and for the length operator 
H where we need the action of the volume operator on non- gauge invariant functions in 
an intermediate step. 



5.3 Examples 

Given a vertex v and a vector of spins j colouring the edges incident at v , let us denote 
by d(j,v) the number of linearly independent vectors \g > compatible with j. This 
number is the dimension of a vector space of vertex contractors associated with v and it 
is given purely algebraically as the number of linearly independent trivial representations 
which appear in the decomposition into irreducibles of the tensor product representation 
7Tj 1 ®..®7Tj n . One can determine this number for each case at hand by repeated application 
of the usual Clebsh-Gordan theorem TCj^nji = ® ■■ ® lx \j-j'\ an d counting all 

the trivial representation of spin that appear. Notice that the sum is direct so that all the 
representations are orthogonal and therefore linearly independent. What one would like 
to have is a closed formula for d(j,v) but that seems to be a complicated combinatorical 
problem. To our knowledge, such a dimension formula is only known for the case of an 
n — fold tensor product of fundamental representations j\ = .. = j n — 1/2 for the groups 
GL(N), U(N), SU(N) |25| . We will not deal with this problem in the present paper. 
First of all, whenever d(j, v ) = 1 then is the only possible eigenvalue (this observation 
was made first in [K| for the case of a trivalent graph in which the vertex contractors 
span only a one-dimensional vector space). More generally, whenever d(j,v) is odd, we 
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know that is an eigenvalue of multiplicity at least one. 

Since M v is a d x d skew matrix its characteristic polynomial has the structure t k (t 2 — 
Xj)..(t 2 — X\d-k)/2) = t h P(d-k)/2(t 2 ) where < k < d for d even and 1 < k < d for d odd. It 
follows from Galois theory that we can determine the spectrum by quadratures exactly in 
general whenever d(j,v) < 9 and in fortunate cases (whenever has multiplicity at least 
d — 8) up to arbitrarily high dimension. In the most general case, however, a numerical 
evaluation is the only possible approach. 

We consider the first non-trivial cases d(j, v) < 9 for an arbitrary valence n of the vertex. 
We may label the spins in such a way that < j± < .. < j n . We may parametrize the 
situation by the no n- negative number k :— ji + .. + j n -i — j n (there is no contractor for 
k < 0). If k is low valued then it is clear from Clebsh-Gordan theory that d(j,v) will be 
low valued. Specifically, if k — then the first n — 1 spins have to add up to maximum 
spin and thus d — 1, the spectrum is only the point irrespective of the value of n. In the 
following we will assume that 2ji > k, otherwise the discussion requires a case division. 
If k — 1 then one sees readily that d — n — 1 : the only irreducible representations 
contained in j 1 <g> j 2 that can possibly add up with j 3 , .., j n _i to j n are j\ + 32, ji + 32 — 1- 
Repeating the argument, the assertion follows. 

Let in general Cfc(n) be the number of representations with weight ji + ..+j n — k contained 
in ji Cg) . . <g)j n for j ra+1 = j\ + • • + jn ~ k to couple to angular momentum where ji <,..,< j n 
and where we assume that j\ — j± < j± + j\ — k, that is k < 2j±. We wish to show that 
Cfc(n) = (n + k — 2, k) where (n, k) is the usual binomial coefficient. Since co(n) = 1 as 
we showed above, the induction is started. Now suppose we know all coefficients up to k 
for each n and start deriving those for k + 1. We can assume that we know them already 
up to n — 1 since Cfe(2) = 1 for all k is trivial to see since we have made the assumption 
2ji > k. Then we have the recursion Ck+i(n— 1) +Cfc(n) = Ck+i(n) which follows from the 
fact that if we have the tensor product j\ <S>..<S>j n -i then there are q(n— 1) representations 
with weight ji + ■■ + j n -i — 1,0 < I < k + 1. But if we now form the tensor product with 
all of them with j n , then the number of times that ji + .. + j n — k arises from those with 
weight j\ + .. + j n -i — 1,0 < I < k, which is given by Ck{n), is equal to the number of 
times that j\ + .. + j n — (k + 1) arises from them. In addition we have Ck+i(n — 1) terms 
coming from the representation j 1 + .. + — (k + 1) itself. It is now readily verified 
that Cfc(ra) solves the recursion. 

Notice that n + 1 is the valence of the vertex. So we get back the result that in the 

tri-valent case we have d = c^(2) = (k, k) = 1 We can now ask for which combinations of 

k, n we have Cfc(ra) < 9. Notice that Cfc(3) = k + 1 : 

Case k = : n arbitrary. 

Case k — 1 : n — 1 < 9 so n < 10. 

Case k = 2 : n(n - l)/2 < 9 so n < 4. 

Case 3 < A; < 8 : n<3. 

Case k > 8 : n < 2. 

Reverse question (only n > 2 makes sense since we need valence at least 3 for the volume 

operator not to vanish trivially) : 

Case n = 2 : k arbitrary. 

Case n = 3: A; + l<9soA;<8. 

Case n = 4 : (k + l)(k + 2)/2 < 9 so k < 2. 

Case 5 < n < 10 : k < 1. 

Case n > 10 : fc = 0. 
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So, for instance we can treat the general 4-valent case under the restriction that j 4 = 
ji + ji + jz — k, < k < 8 for j 3 > j 2 > ji > k/2 in which case d — k + 1. As the formulae 
already get quite involved, we will only discuss the cases k — 1, 2, 3 explicitly. 
Notice first of all that due to gauge-invariance J\ + J2 + J3 + J4IO >= so that we can 
substitute for J 4 . Thus for instance (7124 10 >= —(^121 + ^122 + Q123) 1 > and it is easily 
verified that the operator identity (7122 + 9121 = holds. It follows that 



which simplifies the discussion tremendously since we just need to compute the matrix 
elements of the single operator 5123. We have to compute the k(k + l)/2 matrix elements 



< jl2\qi23\j' 12 >= b'l20'l2 + 1) - j[ 2 U'l2 + !)] H -723 (.723 + 1) < Jl2 | J23 >< jl 2 |j23 > 



for ji2, j'u — ji + 32 — I, I — 0, .., k and the sum extends over j 2 3 = 32 + jz — I, I — 0, ..,k 

which requires the computation of (k + l) 2 6j symbols for j 12 3 — 3 — ji +32 +j3~k = j\. 

The computation is rather tedious. We will display only a few intermediate calculational 

results before giving the eigenvalue A*,. 

Let 2\ 2 := j x + j 2 , J23 := h + h, J = 2(ji + j 2 + ja)- 

Case k — 1. 

The 6j symbols are 



M v = (^) 3 ^[e(l, 2, 3) - e(l, 2, 4) + e(l, 3, 4) - e(2, 3, 4)]g 123 =: aq 12 s 




(5.18) 



and so we have the only non- vanishing matrix element 




(5.19) 



The corresponding eigenvalue of the volume operator is 




(5.20) 



Case k = 2 : 

The 6j symbols are 



< J12IJ23 > = 



j 1 (2j 1 -l)j 3 (2j3-l) 



V'l 2 (2il2 - I)j23(2i23 - 1) 



< 312 - 1 1 J23 > = 



\J j 12 (2j 12 - 2)j 23 (2j 23 - 1) 



< jl2 - 2\j 23 > 



j2 (2j2-l)(J-l)(J-2) 



\ (2j 12 -l)(2jl2-2)j23(2j23-l) 



< .712 1.723 - 1 > 



JiJ 2 (2j 3 -l)(J-l) 

\ Jl2(2jl2-l)j2 3 (2j23-2) 
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< Jl2|j23 
< jl2 ~ l|j23 



2> = 



1 > = 



j 2 (2j 2 -l)(J-l)(J-2) 

\ il2(2ji 2 - l)(2j 23 - l)(2j23 - 2) 

( 2j2 -l)(J-l)-(2j 1 -l)(2j3-l) 
2^'i2(2ji2-2)j 23 (2j 2 3-2) 



< jl2 - 2|j23 - 1 > 

< jl2 - 1 |j23 - 2 > 

-2 > 



(2j 1 -1)(2j 2 -1)j 3 (J-2) 
^(2j 12 -1)(2j 12 -2)j 23 (2j 2 3-2) 



j 1 (2j 2 -l)(2j3-l)(J-2) 
^Ji2(2ji2-2)(2j 2 3-l)(2j 23 -2) 



< jl2 - 2|j23 

and we find then the matrix elements 
< .712 1 9123 Ui2 - 1 > = 



Ji(2j'i-l)j 3 (2j3-l) 



\ (2j 12 - l)(2j 12 - 2)(2j 23 - l)(2j 23 - 2) 



(5.21) 



2j 1 2j 2 (2j 3 -l)(J-l) 1 



x 



< Jl2|gi23|jl2 - 2 > 
< jl2 - 1 1 5l23|jl2 - 2 > 



\ 2(2j 12 -l)(2j 12 -2) 2(2j 23 -i; 
x [(2 3l - l)[6j 3 + 2j 23 - 1] + (2j 2 - 1)[3J + 2j 23 - 7]] 




(2j 1 -1)(2j 2 -1)2j 3 (J-2) 



x 



\ 2(2j 12 )(2j 12 -1) 2(2j 23 -l) 
x [(J - l)[6j 2 + 2j 23 - 1] + (2j 3 - l)[6ji - 2j 23 + 1]] 

(5.22) 

and since the eigenvalues of any 3x3 skew matrix with entries a, b, c are given by 
0, ±i\/a 2 + b 2 + c 2 we find the eigenvalues 0, A 2 where 



^2 = £p\JW\\/< jl2 1 9123^12 - 1 > 2 + < jl2 ~ I|gi23|jl2 ~ 2 > 2 . 



(5.23) 



Case k — 3 i 

We will refrain from displaying the 6j— symbols explicitly in terms of the jj, rather we 
consider them as given through the formulae of the previous section. Now it is easy to 
check that the eigenvalues of a 4 x 4 skew matrix with entries a, b, c, d, e, f are given by i 
times 

t = ±j^Ja 2 + b 2 + c 2 + d 2 + e 2 + f 2 ± |/[^+ f) 2 + (c + d) 2 + (e - f) 2 } x 



x[(a-f) 2 + (c-d) 2 + (e + f) 2 } 

and notice that the argument of the inner square root can be written as (a 2 + b 2 + c 2 + 
d 2 + e 2 + f 2 — 4(a/ + cd — eb) 2 which shows that t is purely real. Substituting the matrix 
elements < ji 2 — / 1 ^123 1 J12 — I' >, </</,< 3 for a, b, c, d, e, f we obtain 



A, 



a \ t 



(5.24) 



The reader should now have a feeling of how to compute the spectrum analytically. The 
formulae derived can serve for analytical estimates as well as the starting point for the im- 
plementation of an appropriate computer code for an algebraic manipulation programme. 
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